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ACCUMULATION RATE OF BOUND STATES OF DIPOLES IN 

GRAPHENE 

SIMONE RADEMACHER AND HEINZ SIEDENTOP 


Abstract. We prove that the bound state energies of the two-dimensional 
massive Dirac operator with dipole type potentials accumulate with exponen¬ 
tial rate at the band edge. In fact we prove a corresponding formula of De 
Martino et al [3]- 


1. Introduction 

Recently the bound state problem for strained graphene with a dipole became of 
interest in the physics literature (De Martino et al [3]). This is described effectively 
by a two-dimensional massive Dirac operator Dfj, with a dipole potential (/>, i.e., 

(1) := a ■ p + as - (j) 

with a := {ai,a 2 ) (the hrst two Pauli matrices), p := (l/i)(i9i, 92 ), and real valued 
potential 

(2) cj) := d + s 
where 



with 5 S is the potential of a pure point dipole at the origin outside the ball of 
radius one around the origin. Without loss of generality, we can - and will from 
now on - pick 0 := (6,0) with 6 > 0, i.e., a multiple of the unit vector along the 
first coordinate axis. The potential s will be the - possibly - singular part of the 
potential that is short range in the sense that — A — s has only finitely many bound 
states. 

It is folklore that the discrete spectrum of would be infinite, if (f) had a 
non-vanishing Coulomb tail. This is also true for its three dimensional analogue. 
However, dimension two and three differ in the case of a dipole potential: whereas 
in three dimensions there are - for small coupling constant - only finitely many 
eigenvalues in the gap [T], ad{D,p) is always inhnite in two dimensions. This has 
been predicted by De Martino et al [3] and proved in [2]. In fact, De Martino et 
al even derived a formula for the accumulation rate of the eigenvalues at the band 
edge. The purpose of this paper is to prove their formula. To formulate our result 
we need some notation: 

Definition 1. We write 

• Ni[A) for the number of eigenvalues of a linear operator A in I G C - 
counting multiplicity, 

• Mh for the - rescaled - Mathieu operator with periodic boundary conditions 
at 0 and 2 tt defined by 

(3) {Mbg){ip) = -g”{p) - bcos{ip)g{p), 
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Bii, for the Bessel operator with imaginary order defined by 


( 4 ) 


{Biuf){z) = -f”{z) - ^f'{z) - ’^fiz). 


Remark 1. The lowest eigenvalue ag of the rescaled Mathieu operator Mf, fulfills 
the transcendental equation (McLachlan [H 3.11.(8)]j 


( 5 ) 


ao= 4 


2 \4J 64 V4/ 576 V4/ 


1 1 Qq 

-L 4 4 


1_1 Qq 

^ 16 4 


1_1 Qq 

^ 36 4 


in particular ao is negative for all b (McLachlan [S] 3.25. diagram]^. 
The above notation allows to formulate our main result. 


Theorem 1. Assume that (j) = d s is real valued with d as in m, and that the 
singular part s is relatively compact with respect to Dq and its negative part s_ even 
fulfills 



s-{x) log(2 + |a;|)da; + 



sl(7rt)| log(i)|]dt < 00 


and also 


s^(x) log(2 + |x|)da; + / (s^)*(7ri)| log(t)|]di < oo. 


Then 


( 6 ) 


lim 

Ey^l 


^(-e,e){D4>) 

|log(l - E)\ 


1 

TT 




Before we embark on the proof we make a few comments: 


Domain: Our condition on the potential (f assures that is self-adjoint on 
: C^) and that the essential spectrum of is M\ (—1,1). 

Electric Potentials: One possible realization of 4> is to think of it as the 
electric potential of some sufficiently smooth and localized charge density 
p, i.e., 



P(>l)dt] 
|(a;,0) -t)| 


with vanishing monopole moment, i.e., p = 0 and to assume that the 
dipole moment r)p(r))dr) of p points in the direction of the first coordi¬ 
nate axis. (See, e.g., Jackson O Chapter 4] for a discussion of multipole 
expansions of potentials. Note that we extend the vector a; G by zero to 
a vector y = {x, 0 ) G K^.) 

Instead of three dimensional densities, we could also allow for densities 
p(y) = P^'^\x)5{xg) that are confined to the electron plane. 

Point Charges: Our hypothesis excludes point charges located directly in 
the graphene sheet, i.e., the plane in which the electrons move. Although 
it is certainly possible to treat a finite number of such singularities with 
subcritical coupling constants, i.e., less than 1/2, this would require an 
analysis of the compactness properties of when restricted to functions 
in a ball containing the singularities of (j>. We refrain from embarking on this 
subtlety, since it is irrelevant for the long range behavior of the potential 
which determines the asymptotic behavior of the eigenvalues. 
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2. Proof of the Main Theorem 


Intuitively the long range behavior determines the asymptotic behavior of the 
eigenvalues. This motivates to prove the following lemma including pure dipole 
potentials. 


Lemma 1. For real a, 6 G R, R € R+, A\= {x € |a;| > 1}. Set 

a + 5a;i/|a;| 
na,b •= —^-i—- 


on Hq{A) or Fl‘^{A), i.e., the operator with Dirichlet respectively Neumann bound¬ 
ary conditions on the boundary of A. Then 


(7) 


^{-oo,E){Ha,b) 

hm —--f—r--- 

E^O |l0g(|S|)| 


271 


trv'(Mf, 


a)_. 


Note that the theorem is only nontrivial, if the sum of the lowest Mathieu eigen¬ 
value —fj, and —a is negative. This, however, is the case for a = 0, since the first 
Mathieu eigenvalue is always negative (McLachlan [5]). 

We know of two ways proving the lemma. The first relies on Dirichlet-Neumann 
bracketing; it adapts an argument that Kirsch and Simon [7] developed for direc¬ 
tion independent pure 1/r^-potentials to the case of dipole potentials (direction 
dependent decay (!)). (See Rademacher [S] for details.) The method presented here 
is somewhat different. As we will see in the proof, the essential ingredient is the 
short range behavior of modified Bessel functions with imaginary order which - we 
feel - is more direct. Moreover it is closer to the derivation by De Martino et al [3]- 
We begin with the proof of the Lemma. 


Proof. We solve = —Xip by separating variables, i.e., we make the ansatz 

ipix) = f{r)g(jp) in spherical coordinates xi = rcos{(p) and X 2 ■= rsin((/3) and 
require periodic boundary conditions for g, i.e., g{ip) = g{ip + 27r). Then g is an 
eigenfunction of the Mathieu operator, i.e., Mbg = —g,g, and / is an eigenfunction 
of the Bessel operator = —A/ which we need to solve with Dirichlet 

respectively Neumann boundary condition at 1 and Dirichlet condition at infinity, 
depending on whether we solve the problem in Hq{A) or H^{A). 

As mentioned above, the Mathieu operator Mb has for all positive b at least one 
negative eigenvalue. Because of the boundary condition at infinity we have that 
the only solutions of (HI are 

(8) fir) = K.^iVXr). 

(Note that the functions li^/ji+Ei'^’) two linearly indepen¬ 

dent solutions of Bessel’s equation. However, li^/p^iVX-) is excluded because of 
its exponential blow up at infinity (Watson [T71 Chapter 7.23, Formula (2)]). 

Next we note that, since the change from Dirichlet to Neumann boundary con¬ 
dition at 1 is a perturbation of rank one for each fixed eigenvalue of the Mathieu 
equation, we merely need to consider the Dirichlet case. Thus we are interested in 
finding the maximal number of nodes that a function Ki^/ji^i'/X-) can have for 
fixed g. a, assuming that = 0 any —A < E, i.e.. 


(9) 


{r>l\K,^iVXr)=0,-X<E} 


If a -|- /r < 0, the operator B has no eigenfunctions and the claim is trivially 
true, i.e., we may assume that a-|-/i > 0. Thus, we need to find the maximal n such 
that 
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where denotes the zeros of using their asymptotic expansion 

(see (IMD). Taking the logarithm and dividing by | log(—i?)| yields 

(11) ^/|iog(-£;)|^ 

as E ^ 0. Since the lower bound can deviate by at most one, this is covered as 
well. □ 


Next we turn to the proof of the theorem: 


Proof. We begin by noting that the operator is a relatively compact perturba¬ 
tion of the free Dirac operator Dq. This implies that has only eigenvalues of 
finite multiplicity in (—1,1), the spectral gap of Dq. Furthermore those eigenvalue 
can only accumulate at —1 or 1. Thus, by the spectral theorem 

(12) fV(_i.i)(D^) = iV(_oo.o)(^^-l)- 

which reduces the problem to study the negative eigenvalues of a relatively compact 
perturbation of the Laplacian, since 

(13) D‘^-1 = -A + (j)'^ + {a ■ p)(t> + (t>{a ■ p) - 2aQ4>- 

The Schwarz inequality - followed by the geometric-arithmetic mean inequality - 
yields for any positive e 

(14) |25R(^/>,e((T-p)(e"V)^)l < ■ 

Thus 

(15) - 1 < (1 + e2)p2 ^ ^ 

and 

(16) e2)p2 + (1 _ g-2)^2 _ 

Note the lower bound (0 is bounded from below for e G (0,1), since both (p and 
0 ^ are relative compact perturbations of p^. 

Both right hand sides of (fT51) and (1161) separate in two independent one compo¬ 
nent operators, since cts is diagonal. We shall focus on the first component. (As 
the proof shows the second component will give the same answer because of the 
symmetry of the pure dipole part.) We write 

(17) (1 ± e^)/i± := (1 ± e^) (p^ + 

for the first components of the right hand side of (Ha and (nni). The task is now 
to estimate N(^_ao,E){h+) from below and A(_oo,£;)(h_) from above as E Z' 0. We 
begin with the lower bound to N(^_ao,o){h+) and write in the spirit of (1^^ 

(18) h± = -A + V± + W±, 

(19) V± : =-2(l±e^)-^d, 

( 20 ) W±: = e-^p^-2(l±e^)-^s. 


(Note that the indices ± at h, V, and W, are just indices motivated by the signs 
in m and HU) and not to the positive part and negative part of an operator as 
elsewhere in the paper.) Thus, by (1^51) 


(21) 


1. N(^-oo,E){h+) 1V(_oo,£:)(-A - (1 

limsup —-j— ,, > hmsup-y- 

Ey'O |l0g(|7^l) E/'O |l0g(|7?|) 
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sinceiV(_oo,£;)(-A + (l-e)e ^W+) < 7V(_^_o)(-A + (1 - e)e ^W+) < oo by (HH) 
and thus vanishes when divided by log(|i?|) as E 0. Finally, we take e to zero 
and get 


( 22 ) 


lim sup 

E/'O 


N{-oo,E){h+) 


iiog(ii?i)i iiog(ii?i)i 

N{-oo,E) ((-A - 2d)|^2({,j,gR2||2,|>l)}^ 


A(_oo,£;)(—A — 2d) 


> lim sup 

E/'O 


|log(|i?|)| 




using © in the last step. Repeating the argument for the second component yields 
the same result. Adding the results for both components gives the claimed upper 
bound. 

We now turn to the upper bound. We use (HU) to estimate the operator from 
below and thus, the number of eigenvalues from above. Again the operator decou¬ 
ples into two one-component operators and we are left with the task to compute 
twice the number of eigenvalues of /i_ below —E as remarked already above. Next 
we use (EZD to estimate from above and note - similarly to the lower bound - that 
the W- part does not contribute. As above we now get 


(23) 


,. A(_oo,_e)(^-) A^(_oo,£;)(—a — 2 d) 

lim sup —-j—< hmsup-^- 

E^o |log(|A|)| eyo |log(|d5|)| 

^ A(_oo,£;) ((-A - 2d)|//2({,rgR2||,r|>i)}) 1, - 

< hmsup ^--j— .|„|'| -= try^Mab. 

Byo log(A) 27r 


where we estimated by the Neumann operator and used © again. Doubling the 
bound because of the two components gives the desired result. □ 
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men). 


Appendix A. Some Auxiliary Formulae about the Asymptotic of 

Bessel Functions 

Dunster [H Formula 2.8] offers the asymptotic formula 

(24) 

= 2 exp(-(n 7 r-())i,)/u) ^1 -f 2(?w^^ (^t.)/u) 0{exp{-A{mT - ^^)/v))'^ 

as n —> 00 for the n — th zero (counting from the right) of the modified Bessel 
function Ku, for fixed imaginary order. Here 

(j)v := arg(r(l -hiu)). 

Appendix B. Bounds on the number of eigenvalues of Schr5dinger 

OPERATORS IN TWO DIMENSIONS 

Since negative potentials in one and two dimensions generate always at least one 
bound state with negative energy, the standard Lieb-Cwickel-Rosenblum bounds 
cannot hold. Following earlier works in the physics literature (see Khuri, Martin, 
and Wu i), Shargorodsky [11] Theorem 4.3] showed that there is a positive con¬ 
stant C such that for all V the number N-{V) := A^(_oo,o)(~A -|- V) of negative 
eigenvalues of — A -|- R is bounded by 

(25) N_{V)<c(^J V-{x)logi2 + \x\)dx +Vl{TTt)\log{t)\]dt^ +1 
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where the subscript V- denotes the negative part of V and V* the spherically 
symmetric rearrangement of V-. 

Appendix C. Known bounds on the number of eigenvalues of sums of 

OPERATORS 

The number of eigenvalues of the sum of two self-adjoint operators A andi? that 
are bounded by below and having infcress(A) = inf dess (^) = 0 can be estimated 
by the sum of the number of eigenvalues of each single operator, i.e., for A < 0 

(26) -^(-oo,£;)(^ + B) < A^(_oo,£:)(^) + N{-oo,e){B)- 

This formula is a consequence of the minimax theorem and was proved by Kirsch 
and Simon [71 Proposition 4] resp. Reed and Simon (TUI p.274, Formula (125)]. In 
particular this formula applies for Schrodinger operators of the form H = —A + 
V + W with potentials V, W in such that inf aess{H) = 0 [71 Proposition 5 (i)]. 
We obtain for e > 0 and E < 0 

(27) iV(_oo,E)(-A + V + W)< + ^F) + + 7 ^^)- 

This estimate leads to the lower bound 

(28) iV(_oo.E)(-A + F + W)> A(_^,E)(-A + (l-e)F)-iV(_^,s)(-A-l^W) 
for the number of eigenvalues of 77 [71 Proposition 5 (ii)]. 
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